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ABSTRACT 

We have compared the halo kinematics traced by globular clusters (GCs) and planetary nebu- 
lae (PNe) for two elliptical galaxies in dense environments (NGC 1399 and NGC 4649) and 
for the merger remnant Cen A. We find differences in the rotational properties of the PNe, 
red GCs, or blue GCs systems in all three galaxies and some indication for substructures. We 
interpret these as signatures of tidal stripping or accretion of small galaxies, showing that low 
luminosity and dwarf galaxies have a higher specific frequency of GCs than PNe. 

Key words: galaxies: elliptical and lenticular, cD; galaxies: haloes; galaxies: individual: 
NGC 1399, NGC 4649, NGC 5128; galaxies: kinematics and dynamics. 



1 INTRODUCTION 

Stellar halos of galaxies offer an important laboratory to understand 
the galaxies' formation process and evolution. The dynamic time 
scales in the halos are large, and the imprint of the formation mech- 
anisms may still be preserved at large radii in the kinematics, in the 
orbital structure, in streams and substructures, or in the chemical 
composition and distribution of stars. 

Due to the rapid dimming of the stellar light at large 
radii, alternate tracers of the underlying stellar population 
are studied. Planetary Nebulae (PNe) and Globular Clus- 
ters (GCs) can be detected in large numbers out to sev- 
eral effective radii, and have b e en used to probe the kine- 



matics (e.a. ICoccato et al. 


2009; 


Herrmann & Ciardullo 
20101: ISchuberth et alj R 


20091; 

mi). 


Teodorescu et al] l20ld: iLee et al.l [ 


dark matter content (e.g. 


Cote et al.l l200ll: iRomanowskv et al. 


2003; de Lorenzi et al.ll200S 


; Napolitano et al. 201 1 ; Richtler et al. 


201 it iDeason et alj 1201 2|) and chemical composition (e.g. 
Forbes et al]|20T it lAlves-Brito et al JI201 1) of the stellar halos in 



nearby galaxies. PNe have been found to be a good tracers of the 
underlying stellar population, as their sp atial distribution and kine- 
matics well agree with those of the stars dCoccato et alj2009t) . This 
is still not so clear for GCs: they have a bim odal color distribution 
(see the review by iBrodie & Strader 2006 and references therein) 
where red GCs (metal rich) have a radial density distribution simi- 
lar to that of the stars and PNe, while blue GCs (metal poor) have 
a more extended distribution. 

Recently, the distribution and kinematics of PNe and GCs 
in galaxy halos have been exploited also to detect substructures, 
which are interpreted as signatures of minor mergers and accre- 
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tion events that contri bute to the build up of the stellar halos 
in ea r ly-type galaxies {Shih & Mended l20ld : IWoodlev & Harris! 
l201lt Romanowskv e t alj|2012f) or to the disk heating in spirals 
dHerrmann et alj2009bl) . 

So far, a direct and detailed comparison between the two- 
dimensional velocity fields inferred by GCs or PNe has been done 
only in few cases. The commonly used approach was either to com- 
bine the two datasets in the analysis, or use the most suited dataset 
to study a particular science case. Little attention was put in the 
information potentially contained in the "discrepancies" between 
the PNe and GCs velocity fields. We therefore compared the halo 
kinematics independently traced by GCs and PNe in a sample of 
galaxies. In this paper, we describe the properties of 3 galaxies in 
our sample (namely NGC 1399, NGC 4649 and NGC 5128), whose 
two-dimensional velocity fields have kinematic quantities that are 
different for PNe and GCs. The complete analysis of the entire sam- 
ple will be presented in a forthcoming paper. Here we concentrate 
on the average kinematic properties and we do not discuss the radial 
distribution of individual tracers. 



2 HALO KINEMATICS 

We reconstruct the mean two-dimensional velocity fields < V > 
for the PNe, red-GCs and blue-GCs systems using an adaptive 
Gaussian kernel smooth i ng te chnique, which was defined and 
tested in ICoccato et ID J2009). 3a velocity outlier s are itera- 
tively removed through the friendless algorithm as in lMerrett et al.l 
(2006). We do not assume any particular symmetry of the data to 
not bias the analysis towards any preferential configuration or sym- 
metry. The PNe and GCs detected in NGC 1399 and NGC 4649 
contain objects associated to the companion galaxies (NGC 1404 
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Figure 1. Membership allocation of PNe and GCs between the main galaxies (NGC 1399 and NGC 4697) and their companion (NGC 1404 and NGC 4647). 
The first figure in each row shows the radial velocity of PNe and GCs as a function of their position along the line connecting the main galaxy to the companion. 
The second and third figures in each row show the location of the PNe and GCs on the sky. Objects marked with red have more than 60% probability to belong 
to the companion galaxy, and therefore are removed from the analysis. Green triangles mark the position of the galaxies and the companions. 



and NGC 4647, respectively). Therefore, we have to remove the 
contaminants before recovering the velocity fields. 

T o do that, we use the procedure described in iMcNeil et al.l 
which assigns to each object (PNe or GCs) a probability 
to belong to the galaxy or to the companion on the basis of the 
galaxies' surface brightnesses at the objects' positions and of the 
line of sight velocity distributions. We removed from the analysis 
all objects that have more than 60% of probability to belong to the 
companion galaxy (see Figure [T}. This is different from what was 
typically used as criterion in the past (e.g. objects with more than 
95% probability to belong to the main galaxy were kept) because 
we want to include in the analysis also those objects that may have 
been stripped from the companion and thus trace an on-going in- 
teraction. The reconstructed mean two-dimensional velocity fields 
of PNe, red and blue GCs are shown in Figure[2] 

We measure the amount and direction of rotation as a function 
of radius for each of the components (PNe, red and blue GCs) in 
each galaxy by fitting the following function for several radial bins: 

< V > (0; R) = V sys {R) + KotOR) ■ cos \<j> - PA kin (R)] (1) 

R is the distance of the object (PNe or GCs) from the galaxy 
center, cj> is its position angle on the sky, < V > is the mean ve- 
locity of the objects at a given position the sky. The fit variables 
for each radial bin are: Kot (amplitude of rotation), V sys (mean 
velocity), and PA^n (kinematic position angle). V aya is found to 
be consistent with the galaxy systemic velocity at all radii. Angles 
are measured from North towards East. Errors on the measured fit 
variables V sys , Kot, and PAkm, are computed by means of Monte 
Carlo simulations. We create 500 mock catalogs with simulated 
velocities V{,qs for each object drawn from a Gaussian distribu- 
tion with mean VLos and standard deviation that includes the lo- 



cal velocity dispersion and measurement errors. We then reprocess 
the simulated catalogs and get 500 simulated fit parameters, whose 
standard deviations define the errors of our fit variables. 

In the following, we describe the results obtained for the 3 
galaxies individually. 

NG C 1399. PNe data are f rom lMcNeil etal] ( l2010l) , GCs data 
arc from lSchuberth et al.l d2O10b . 

The PNe system is dominated by random motions (< a >~ 
220 km s _1 ), but it shows also small rotation along two opposite 
directions. Innermost R < 200" PNe rotate slowly (Vrot = 31 ± 
20 km s" 1 ) along PA kin = 231° ± 50°, co nsistent with the th e 
galaxy photometric major axis (PA = 290° , iDirsch et al]|2003l) . 
On the contrary, outermost R > 300" PNe rotate slightly faster 
(V^t = 47 ± 27 km s _1 ) and along a nearly opposite direction 
(Pylkin = 46° ± 47°). The PN e two-dimensional v elocity field is 
in agreement with that shown in IMcNeil et al]d2010b - 

We divide t he GC s in NGC 1399 in red and blue GCs as in 
ISchuberth etall j2010l) . They find a rotation of V = 110 ± 53 
km s _1 for the outermost (4' < R < 8') blue-GCs. From 
Equation [TJ we find consistent values in the same radial range 
(V^t = 79 ± 29 km s _1 ). Nevertheless, by looking at Fig. [2] 
both the red and blue GCs do not display a well defined direction 
of rotation: their mean velocity fields are rather characterized by 
substructures. The two more evident are located on the South and 
North East sides of NGC 1399. In Figure [3] we show the line of 
sight velocity distributions in these regions. In the S substructure, 
the PNe and GCs have median velo cities that differ a t ~ 1.2a level. 
From a Kolmogorov-Smirnov test jPress et ail2002l) . the probabil- 
ity that PNe and red (blue) GCs are drawn from the same distribu- 
tion is 1 ± 1% (8 ± 5%); this probability is higher when comparing 
the red and blue GCs distributions (69 ± 22%). In the NE substruc- 
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Figure 2. Mean velocity fields for NGC 1399 (top panels), NGC 4649 (middle panels) and NGC 5128 (bottom panels). Each panel shows the galaxy image, 
where the positions of the PNe and GCs are marked with crosses. The colors indicate the values of the mean velocity field computed at that position. North is 
up, East is left. Coordinates are in arcseconds. The dashed lines define regions where kinematic substructures in NGC 1399 are identified (see text for details). 



ture, the median velocities of red and blue GCs differ at ~ 1.2a 
level (no PNe are detected in this region), and the probability that 
they are drawn from the same distribution is 3 ± 2%. Note that 
Fig.fJJshows Vlos, which are independent of the kernel smoothing 
technique used to recover the mean velocity fields. 

NGC 4 649. PNe data are from lTeodorescu et alj l l201ll) . GCs 
data are from lLee et al, (2008). 

The PNe system is dominated by dispersion (< a >~ 250 
km s -1 ) and it shows two main kinematic substructures. PNe in the 
inner R < 60" rotate with an amplitude V ro t = 77 ± 22 km s~ 
along PAun = 302 ° ± 20° , consistent with the photometric major 
axis (PA = 285°, |de Vaucouleurs et a PNe in the outer 

R > 250" rotate around an orthogonal direction (PA^in = 205° ± 
18°) with the same amplitude (Kot = 75 ± 29 km s" 1 ). The 
two-di mensional velocity field is similar to that shown bv lDas et al.l 
J201 l|), obtained using different criteria for contaminants removal. 

We divide the G Cs system of NGC 4649 in red and blue GCs 
as in Lee et al. ( 2008). The velocity field of the red GCs is not easily 
described by Eqn.Q] therefore it is not possible to define an ampli- 



tude and direction of rotation. The blue GCs system shows a more 
regular velocity field. In the inner R < 60" the blue GCs rotate 
along PA kin = 292° ± 31° with an amplitude of Vrot = 33 ± 17 
km s _1 , slower than the PNe in the same radial range but along 
similar direction. The PAkin and V TO t increase gradually with the 
radius. At R > 300" we measure rotation V I0 t = 82 ± 20 km s _1 
along P^kin = 342° ± 88°, nearly opposite to that of the outer- 
most PNe. 

The stars in NGC 4649 within ~ 25" rotate along a direc- 
tion very similar to those of the inner PNe and blu e GCs (PA t = 
271.5 ° ±3.8°) with an amplitude of ~ 94 km s" 1 dKrainovic et al.l 
1201 ll) . The stellar rotation agrees with that of PNe, but it is larger 
than that of the blue GCs. This difference is probably due to the 
spatial incompleteness in the center that overestimates the kernel 
amplitude adopted in this region, causing a dilution effect when 
computing the local mean velocity field. 

The dependence on radius of Vrot and PAkin for stars, PNe, 
and blue GCs in NGC 4649 are shown in Figure [4] 
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Figure 3. Velocity distribution of PNe (green), red GCs (red) and blue GCs 
(blue) in the South (left panel) and North East (right panel) substructures of 
NGC 1399. Arrows mark the median radial velocities of the systems with 
the same colors. Velocities are in the heliocentric reference system. 



NG C 5128. PNe data are from IPeng et all d2004ij) . GCs data 
are from lWoodlev et all J2010al lbh- 

The PNe and GCs kinema tics i n NGC 5128 were previousl y 
compared by dPeng et alj2004bl) and lWoodlev et aljj200l[2010bh . 
Here, we re-analyze these data sets uniformly with respect to those 
of the other galaxies in our sample, find i ng con sistent results with 
their studies. Following IWoodlev et alj feOlObl) , we use the GCs 
metallicity [Fe/H] = — 1 as threshold to separate red (metal rich) 
and blue (metal poor) GCs. 

The PNe velocity field shows disk-like rotation, with a twist 
in the central 5 k pc due to the triaxiality of the system (see also 
IPeng et al]2 004a). The small asymmetry between the directions of 
the approaching and receding side of the galaxy are probably due to 
the non uniform sampling of PNe in the outer regions. From Eq.Q] 
we find that PNe reach a maximum rotation Kot = STilSkms^ 1 
around 20', and then a decline to V TO t = 56 ± 13 km s _1 in the 
halo (consistent with lWoodlev et alj2007l) . The velocity dispersion 
is ~ 140 km s _1 in the center and then it declines with radius 
down to ~ 80 km s _1 . 

GCs do not show a clear disk-like rotation as PNe, especially 
the blue GCs. Central velocity dispersion of GCs is consistent with 
that of PNe, but it does not decline with radius. This support a pic- 
ture where the GCs are more related to a spheroidal component, 
while the PNe are more related to a disk-like component. 



3 DISCUSSION AND CONCLUSIONS 

We showed in Sect. |2] that the PNe and GCs systems in NGC 
1399, NGC 4649 and NGC 5128 have different kinematic prop- 
erties. NGC 1399 and NGC 4649 are located in dense environ- 
ments and have a luminous companion, respectively NGC 1404 and 
NGC 4647, while NGC 5128 is considere d a merger remnant (e.g. 
IPeng et alj|2004bl ; IWoodlev et alj|2010al and references therein). 
Therefore, the observed discrepancy between PNe and GCs kine- 
matics in these galaxies could be related to the environment or to 
the merging history. We now discuss some scenarios: 

Tidal stripping. NGC 1399 and NGC 4649 could have stripped 
material from their less massive companions. The GCs system 
(the blue GCs in particular) is ge nerally more extended than the 
underlying stellar population (e.g. lHarris|[T99ll) : this implies that 
in some cases only the outer GCs will be stripped by the tidal 
interaction, and no stars and PNe. If this is the case, the sig- 
nature of the tidal stripping will be visible preferentially in the 
blue GCs system, which will therefore trace different kinemat- 
ics than the PNe and stars. iForbes et al] dl997l) explored this sce- 
nario for NGC 1399. This galaxy has an high specific frequency 




100 



200 300 
R [arcsec] 



400 



Figure 4. Rotation amplitude (V TO t) and direction (PA^^) of the stars 
(black star), PNe (green), and blue GCs (blue) systems in NGC 4649, as 
obtained from EquationfT] The dotted lines indicate the photometric major 
axis position angle. 



(i.e. number of GCs per unit galaxy lum inosity, normalized at 
M v = 15, lHarris & van den Berghlll98ll) , and GCs are prefer- 
entially distributed in the outer galaxy regions. This suggests that 
some GCs in NGC 1399 may come from NGC 1404, which has 
a low specific frequency, consistent with having lost many of its 
GCs. Moreover, NGC 1404 has an observed luminosity which is 
consistent with its stellar velocity dispersion, suggesting that only 
a small fraction of its stars were removed in the stripping process. 

Within the GCs system, the hypothesis that blue GCs can trace 
the interaction between the two galaxies more clearly than the red 
GCs system is supported by the following considerations: i) the ma- 
jority of GCs with redshifted velocities in NGC 4649 are blue, and 
are located on the NGC 4647 side (Figs.Q]and[2]l; and ii) the blue 
GCs in the South substructure of NGC 1399 (close to its compan- 
ion) have a broader line of sight velocity distribution than the red 
GCs in the same region. 

Kinematic disturbances from a companion. Kinematic distur- 
bances induced by the gravitational interaction with the companion 
can generate asymmetries in the stellar and ionized- gas velocity 
fields as observed in interacting bina ry galaxies (e.g. Borne et al.l 



[1991) and galaxies in clusters (e.g. iRubin et al.H l999). However, 
stars, PNe and GCs would feel the same gravitational influence 
from the companion, and therefore they would trace the same ve- 
locity fields. Thus, we exclude this scenario. 

Multiple mergers and accretion episodes. Massive galaxies 
in dense environments are subje ct to many merger and accre- 
tion events of different kinds (e.g. |Pe Lucia & Blaizoll2007l) . Al- 
though NGC 5128 does not live in an environment as dense as 
those of NGC 1399 and NGC 4649, the presence of substruc- 
tures ( Woo dlev & Harrisl201 lljMouhcine et al.l201 1) suggests that 
NGC 5128 experienced significant accretion episodes in the past, 
perhaps from late infall of tidal tail material after the merger. If 
the fractions of accreted GCs and stars/PNe differ, the kinematic 
signature would be more visible in one system than in the other. 
It is therefore interesting to compare the expected number of ac- 
creted GCs and PNe for various galaxies to see whether there is a 
particular class of galaxies that would leave a kinematic signature 
preferentially in GCs or PNe. In Figure[5]we compare the expected 
number of accreted GCs and PNe as a function of the total magni- 
tude and morphological type of the accreted galaxy. Despite of the 
scatter, the low statistics for dwarf galaxies, and the different obser- 
vational constraints for PNe and GCs, there is indication that low 
luminosity and dwarf galaxies contribute more GCs than PNe. On 
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Figure 5. Expected number of sources (GCS or PNe) as a function of 
galax y absolute m agnitude (left) and morphology (right). GCS data are 
from Brodie & Strader 12006) and references therein. PNe data are from 
IZiilstra et alj fe006t) . iBuzzoni et alj fe006t) . ICoccato et all 120091) . and ref- 
erences therein. The expected numbe r of GCs and PNe are com puted from 
i) the specific frequencies, Sn lHarris & van den Bergh|[T98ll) and «i,b 
(Jacoby 1980), for GCs and PNe respectively, and ii) the total luminosities 
in the V and B bands, for GCs and PNe respectively . 



the contrary, the expected contributions from bright spiral galaxies 
are more similar. Therefore, when a massive galaxy accretes many 
low luminosity galaxies and dwarfs, it will preferentially accrete 
GCs. Thus, the signature of the merging events will be imprinted 
in the GCs kinematics system more strongly, and differences in the 
velocity fields of GCs and PNe are to be expected. 

To summarize, there is good indication that the velocity fields 
of PNe and GCs are different; to explore this better will require 
large homogeneous samples of both tracers. Such differences can 
be interpreted as a signature of accretion of dwarf galaxies or 
tidal stripping of GCs, which are more loosly bound than PNe. 
It would be interesting to compare the PNe and GCs kinemat- 
ics in other massive galaxi es in dense envi ronments, like NGC 



4889 JCoccato et al. 



Arnaboldi et al _ _ 
1316 dRichtler et alj 



]2010j | bT), NGC 3311 JCoccato et alj l201ll ; 
20121). M87 TRomanowskv et alj|2012l) . and NGC 



where the evidence for accretion of 



smaller galaxies comes from independent arguments and tech- 
niques. 
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